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We analyze the proposal of achieving a Mott state of Laughlin wave functions in an optical lattice 
[M. Popp et al., Phys. Rev. A 70, 053612 (2004)] and study the consequences of considering the 
anharmonic corrections to each single site potential expansion that were not taken into account until 
now. Our result is that, although the anharmonic correction reduces the maximum frequency at 
which the system can rotate before the atoms escape from each site (centrifugal limit), the Laughlin 
state can still be achieved for a small number of particles and a realistic value of the laser intensity. 
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I. INTRODUCTION 

The fractional quantum Hall effect (FQHE) is one 
of the most studied phenomena in Condensed Matter 
Physics ill . Despite the fact that a complete understand- 
ing of it is still missing, it is commonly believed that the 
interactions between the particles are essentially respon- 
sible for the strange states of matter that the 2D elec- 
tron gas shows at some particular values of the trans- 
verse magnetic field. In this respect, in 1983, Laughlin 
proposed an Ansatz for the wave function of the ground 
state of the system [2|. This wave function, defined by 

* m (zi, . . . , ZN ) ~JJ(zi - Zi )- e -sf =1 i^r/2 1 (1) 

where Zj = Xj + , j = 1, . . . ,N stands for the position 
of the j-th particle, and m is an integer number, describes 
the fractional quantum Hall state at a filling fraction v = 
1/m. However, this state has only been proven to be 
an exact eigenstate of very specific Hamiltonians Q and 
for some specific values of the filling fraction; it contains 
the relevant properties that the ground state of the real 
system must have. 

One of the most important features of some of the 
fractional quantum Hall states is that they are states of 
matter whose quasiparticle excitations are neither bosons 
nor fermions, but particles known as non-Abelian anyons, 
meaning that they obey non-Abelian braiding statistics. 
These new phases of matter define a new kind of order in 
Nature, a topological order Q • Such systems have become 
very interesting from the quantum computation perspec- 
tive. Quantum information could be stored in states with 
multiple quasiparticles, which would have a topological 
degeneracy. The unitary gate operations would be simply 
carried out by braiding quasiparticles, and then measur- 
ing the multi-quasiparticle states. In this respect, let us 
also mention that several spin systems with topological 
order has been proposed recently as candidates for robust 
quantum computing 0, 0] . It has also been shown that 
these proposals could be implemented by means of cold 
atoms 0,111 ■ 

Despite the great possibilities of the FQH states, nei- 



ther the Laughlin wave function nor its anyonic excita- 
tions have been observed directly in an experiment so far. 
Nevertheless, recent experimental advances in the field of 
ultra-cold atomic gases suggest that they could be good 
systems to simulate many Condensed Matter phenom- 
ena, and, in particular, the FQHE. From the theoretical 
point of view, it has been shown that the FQH can be re- 
alized by rotating a bosonic cloud in a harmonic trap (see 
Ref. 0, E3|). The rotation plays the role of the magnetic 
field for the neutral atoms. Thus, in the fast rotation 
regime, the atoms live in the lowest Landau level and, if a 
repulsive interaction is introduced, they form the Laugh- 
lin wave function. In such system, the Laughlin state is 
a stable ground state, however, in practice, due to the 
weak interaction between the particles, the gap is too 
small, and it is not possible to achieve it experimentally. 

An idea to avoid this problem is to use optical lattices 
[ill [l2l ]. In such systems, the interaction energies are 
larger, since the atoms are confined in a smaller volume. 
This yields a larger energy gap and, therefore, opens the 
possibility of achieving FQH states experimentally. 

In the present work, we improve the existing proposal 
of achieving the Laughlin state in an optical lattice pre- 
sented in Ref. [l3[ by studying the consequences of con- 
sidering not only the harmonic approximation to the sin- 
gle site potential expansion, but also the anharmonic cor- 
rections. Our work is organized as follows: first, we ana- 
lyze the one body Hamiltonian of our system. The solu- 
tion of the 2D harmonic oscillator is reviewed and a more 
realistic model for the single site potential of a triangu- 
lar optical lattice is discussed. Then, we study the many 
particle problem by introducing a repulsive contact in- 
teraction between the particles. We compute, by means 
of the exact diagonalization technique, the fidelity of the 
ground state of the system with the Laughlin wave func- 
tion for a wide range of the experimental parameters. 
Finally, we discuss which are the required experimen- 
tal conditions and the procedure to obtain the Laughlin 
state. 
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II. ONE BODY HAMILTONIAN 

A. Harmonic case 

Let us consider one atom confined in a harmonic po- 
tential which rotates in the x — y plane at a frequency O. 
We will assume that the confinement in the z direction is 
sufficiently strong so that we can ignore the excitations 
in that direction. The Hamiltonian associated with this 
system is 
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H '° = 2M & + ® + 2 W(f2 + y 2) - nL > (2) 

where M is the mass of the particle, p x and p y are the 
canonical momenta associated with the position coordi- 
nates x and y, uj is the frequency of the trap and L is the 
angular momentum operator. 

We can easily diagonalize this Hamiltonian by defining 
the creation and annihilation operators of some circular 
rotation modes a± — (a x ± ia y ), where the 
standard annihilation operators of the ID harmonic os- 
cillator in the directions x, y. In terms of the number 
operators n± = a ± d± corresponding to these circular cre- 
ation and annihilation operators, the angular momentum 
can be written as L = ft (n+ — ri_ ) and the Hamiltonian 
reads 



huj(u) — fl)h^ 
Its spectrum is, therefore, 



+ H(u> + f2)n_ + hjj . 



(3) 



n+ ,n. 



= h(oj - fi)n + + H(uj + Q)n_ + Hoj , (4) 



where n± are the integer eigenvalues of the number op- 
erators h±. 

In the fast rotation regime (il — > oj), the family of 
states 



0) 



|0) 

/ 977, \ / 



(5) 



where |0) is the state annihilated by both a+ and a_, 
form the subspace of lowest energy, usually known as 
the lowest Landau Level (LLL). We can write the wave 
function of these states 
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(6) 



where \z) — \x, y) are the eigenstates of x and y, z = x+ ^ v 

is a complex variable and I = \J is a characteristic 
length of the system. 

Notice that if the rotation frequency is high but smaller 
than the trap frequency, the lowest energy subspace is 
only formed by states with m < In this limit, 

we can find an effective Hamiltonian of our system by 
projecting the original Hamiltonian onto the LLL, 



H = h{u - fl)L . 



(7) 



Indeed, although Hamiltonians of Eqs. ([3]) and ([7]) 
are different, their projections onto the LLL coincide, 
PH' Q P^ = PH a P\ where P = J2 m \ m )( m \ is 

the pro- 
jector onto the LLL. 

In Ref . [13[ , the formation of fractional quantum Hall 
states in rotating optical lattices is studied under the ap- 
proximation of a single particle Hamiltonian as the pre- 
sented in Eq. ([7]). Nevertheless, the potential of a site 
of an optical lattice is not infinite as the parabolic one. 
Next, we model in a more realistic way this single site 
potential and study how it will affect the formation of 
fractional quantum Hall states. 



B. Quartic correction 

1. From the lattice potential to our model 

We will consider a system of bosonic atoms loaded 
in a 2D triangular optical lattice The lattice is 

2D due to the confinement in the x — y plane cre- 
ated by two counter propagating lasers in the z direc- 
tion. In each of this planes a triangular lattice is re- 
alized by means of 3 lasers pointing at the directions 
ki = (sin6>,0,cos6>), k 2 = (- sin0/2, %/3 sin6>/2, cos 6>) 
and &3 = (— sin 0/2,— y3 sin 9/2, cos 0), where 9 is the 
angle between these directions and the z axis. 

In the region where the atoms are, the electric field 
created by each of these lasers can be modeled as a plane 
wave, Ej(x,t) = Aee l( - k j' s ~ ULt \ where kj — ^-kj for j 
1,2, .'! are the wave vectors of the lasers and A is the 
wavelength of the light used. Notice that all the lasers 
share the same frequency ujl and polarization vector e. 
The atoms are subject to the superposition of the electric 
field Ei(x,t) created by each of the lasers, 



E(x, t) = Y^ t) = E (x)t 



(8) 



After averaging over time, the effective potential that 
atoms feel is, therefore, the correction to the energy of 
its internal state due to the AC-Stark effect. This energy 
shift is proportional to the square of the amplitude of the 
electric field |i?(af)| 2 and, in our particular case, it can be 
written as 



(9) 



where the sum runs over the 3 different couples of 
ij = 12,23,31, and Uo is the intensity of the laser. 
The minima of this potential form a Bravais triangu- 
lar lattice with basis vectors a\ = fj-^-g i y ~\-> an< ^ 
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f (i,o). 



Let us assume that the intensity of the laser U S> E R , 
where Er — h 2 k 2 /2M is the recoil energy, with fc = 
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2ir/\. In this limit, tunneling of atoms between different 
sites is forbidden, and the lattice can be treated as a sys- 
tem of independent wells. Thus, the ground state of the 
system is a product state of the state of each site (Mott 
phase) and we can study the whole system by studying 
each site of the lattice independently of the others. 

In order to describe the potential of a single site, the 
previous lattice potential in Eq. ^ is expanded around 
the equilibrium position x = 0, 

U{x) ~ ~Mu,V + y 2 ) - ^U (x 2 + y 2 ) 2 + 0(x e ) . 

( 10 ) 

We get a harmonic oscillator term with frequency u> — 
2y/UoEn/h plus a quartic perturbation and higher or- 
der corrections. Notice that both terms of the expansion 
respect the rotational symmetry. This is actually the rea- 
son why we have considered a triangular lattice instead 
of the simpler square one. In the square lattice, the circu- 
lar symmetry is broken in fourth order of the expansion, 
while in the triangular one, this does not happen until 
the sixth order. From now on, it will be considered that 
our single site potential is adequately described by only 
these harmonic and quartic terms of the expansion of the 
potential. This approximation is reasonable since it has 
already been assumed Uq ^> En- 

Furthermore, by introducing some phase modulators 
into the lasers that form the lattice, it is possible to cre- 
ate time averaged potentials that generate an effective 
rotation of the single lattice sites. If this rotation 51 is 
close to u, the system is in the fast rotation regime and 
we can obtain the low energy Hamiltonian proceeding as 
before. Written in units of fku, it becomes 

where 7 = jj^ U °^f — S3\/ Th IS a P ei "turbation param- 
eter. Notice that the perturbation parameter, 7, only 
depends on the intensity of the laser in units of the en- 
ergy recoil. 

In summary, we have described the whole sophisticated 
optical lattice by a set of independent wells modeled by 
the simple Hamiltonian presented in Eq. (TTTT) . This effec- 
tive Hamiltonian is the same from Eq. (J7]) plus a quartic 
correction term which will be responsible for all the new 
physics that are going to be studied next. 



2. Exact solution 

Next, we find the lowest energy eigenstates of the 
Hamiltonian presented in Eq. (jTTJ) . First of all, we realize 
that both H and Hq are rotationally invariant, therefore, 
[if, H ] = [if, L] = 0. Thus, if, if and L must have a 
common eigenbasis, and this can only be {|m)}. We de- 
termine the eigenvalues of if by computing the expected 



values 

E m = (m\H\m) = fl - m-j{m + l)(m + 2) , (12) 

in units of hu. 

Let us point out that the commutator between the 
Hamiltonians of the harmonic and anharmonic systems 
before being projected onto LLL, [H',H' ], is not zero. 
Therefore, their lowest energy eigenstates only coincide 
in the fast rotation regime. 

It is important, then, to establish which is the fast ro- 
tation regime for the Hamiltonian (|11[) or, in other words, 
under which conditions our low energy description of the 
system is correct. To see this, we can compute, using 
perturbation theory, the first correction to the energy of 
the second Landau level (NL) states 

£^ i = (m + 2)-7(m + 7)(m + 2). (13) 

We realize that the LLL is a good description of the low- 
est energy states, i. e. E m <C E^ L , if 37(771+ 1) <C 1. In 
this regime, the projection of the Hamiltonian onto the 
LLL is a good effective description of our system. From 
now on, this condition will be always assumed. 



C. Maximum rotation frequency 

In the harmonic case, we realize that if the rotating 
frequency f2 exceeds the frequency of the trap uj (see 
Eq. ([?]))■ the lower energy states have an infinite angu- 
lar momentum, i. e. all the particles of the system are 
expelled from the trap. This maximum frequency is usu- 
ally known as the centrifugal limit. What we are going 
to study now is, therefore, how the introduction of the 
quartic perturbation will affect the centrifugal limit of 
the system. In order to answer this question we realize 
two different analysis that give quite similar results. 

The first argument consists of a semi-classical interpre- 
tation of the dynamics of the particle. In the anharmonic 
system, we have a competition between the attractive 
force of the harmonic trap and the repulsive forces cor- 
responding to the quartic correction and the fictitious 
centrifugal force. In the central region, the harmonic 
trap dominates and the particles are confined, while if a 
particle was in the exterior region, it would be expelled. 
The limit radius between these two regions corresponds 
to the maximum of the effective potential that includes 
both the trap and centrifugal terms. Furthermore, from 
the semi-classical point of view, each stationary state \m) 
follows a circular trajectory of radius r ~ \pm around the 
origin. We expect then that those states whose associ- 
ated radius is less than the limit radius are bound states. 
Thus, we find that the maximum rotation frequency, or 
the centrifugal limit, depends on the maximum angular 
momentum that we want to keep in the trap. This is a 
remarkable difference with respect to the harmonic case, 
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FIG. 1: Spectra of the single-particle Hamiltonian shown in 
Eq. that models the dynamics of a particle in a rotat- 

ing single site of a triangular optical lattice. This plot has 
been realized for 7 = 0.01 and SI = 0.8w. Notice that for 
these particular values of the frequency rotation and the laser 
intensity, particles with an angular momentum larger than 
rriL — 8 would be expelled. 



where the centrifugal limit was the same for any angu- 
lar momentum state. In particular and according to this 
semi-classical criterion, the centrifugal limit for a particle 
with an angular momentum equal or less than uil is 



Cl L = Uy/l - AjniL ~ uj(1 - 2 7 m L ) 



(14) 



Another approach is to think in terms of the slope of 
the spectra E m presented in Eq. (fT2j) and plotted in Fig. 
[TJ We realize that the single particle energy E m begins 
growing with m from m = 0, reaches a maximum at 



+ 



(where the brackets [ ] represent the in- 
and becomes monotonically decreas- 



27 1 2 
teger part function) 
ing from this point. We interpret this in the same way as 
the harmonic case, that is to say, those states \m) that are 
in an energy interval with a negative slope are unstable. 
According to this interpretation, the angular momentum 



2 7 



in order 



of a particle must be less than 

to be trapped. This implies that the centrifugal limit for 
a trapped particle with angular momentum is 



Q, L = w(l - 2 7 m L - 4 7 ) . 



(15) 



Although the previous approaches are different, notice 
that they predict practically the same centrifugal limit. 
In our simulations presented in Sec. HV| we have taken the 
condition given by Eq. (|15j) . since it is the most restrictive 
one. 

At this point, we have completely solved the one par- 
ticle problem. Let us note that whereas in the harmonic 
case the maximum rotation frequency coincides always 
with the trap frequency, in the system with the quartic 
perturbation, this centrifugal limit depends on which is 



the maximum angular momentum of a particle that we 
want to keep trapped. This difference will be crucial to 
understand why it will be more difficult to obtain the 
Laughlin state in the anharmonic case than in the har- 
monic one. 



III. MANY PARTICLE PROBLEM 

The Laughlin state is a strongly correlated state, there- 
fore, interaction between atoms will play an essential role 
in its experimental realization. 

Let us consider, then, a set of N bosons trapped in a 
well as the one described previously, and interacting by 
means of a repulsive contact potential, 



Vif-r 1 ) = gS{r-r > ). 



(16) 



where the parameter g accounts for the strength of the 
interaction, and is related to the s-wave scattering length, 
a s , and to the localization length in the z direction, £ z , 

We are interested in knowing the ground state of the 
system for a wide range of g and Uq/ Er in order to see 
under which conditions the Laughlin state could be re- 
alized experimentally. The solution of this problem is 
trivial in the extreme cases g — > and g — > 00. 

When g = 0, the atoms do not interact and the ground 
state of the system is merely a product state with all the 
atoms with angular momentum zero. 

If g — > 00, the atoms will find a configuration in which 
they are always in a different position with the mini- 
mum total angular momentum possible. This configura- 
tion that minimizes the interaction energy at the expense 
of the angular momentum of the particles is precisely the 
Laughlin state. As we can see in Eq. ([T]), the interaction 
energy in the Laughlin wave function is strictly zero since 
the probability that two atoms are in the same position 
is null. 

For a finite g, the problem of finding the ground state 
of the system has to be addressed numerically. In partic- 
ular, we will solve it by means of exact diagonalization. 
With this aim, we can write the Hamiltonian of this sys- 
tem in second quantization form, 



H = Y, ((1 - 37 - -7h 2 i+l- 2 7 ) 



(17) 



where 



and bj are the creation and annihilation oper- 



ators of the mode j which corresponds to the LLL state 
with angular momentum j presented in Eq. ^ , hj = bj bj 
counts the number of particles with angular momentum 
j, and Vijki are the coefficients of the interaction in this 
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basis defined by 

V m i m,2 1TI3 rri4 — 

(mi m2|F|m3 7714) (18) 
(mi + m 2 )! 
~ 2tt % /^i!^2!m 3 !m 4 ! 2 m ^+ m ^ 

The cylindrical symmetry of the Hamiltonian allows the 
diagonalization to be performed in different sub-spaces of 
well defined total z-component of L = X^iLi TO *- 

Thus, given the parameters f2, Uq/Er and g and a sub- 
space of total angular momentum L, we construct the 
multi-particle basis of N particles compatible with L, 
calculate the matrix-elements of the Hamiltonian in this 
subspace by means of Eq. (fTT|) . and perform its diago- 
nalization. Once the ground state for each subspace L 
is determined, we find the ground state of the system by 
selecting the state with lowest energy. 

An important issue that we have to take into account 
when we perform the simulation is to be careful with 
the centrifugal limit. As it has been shown in Sec. Ill CI 
if we want that the particle with the maximum angular 
momentum does not escape from the trap, we have to 
keep the rotation below the centrifugal limit given by 
Eq. (I15p . In our case, in which we want to drive the 
system into the Laughlin state, this centrifugal limit is 
determined by the maximum angular momentum of a 
single particle in the Laughlin state, that is N(N — l)/2. 
Thus, the maximum rotation frequency is given by Eq. 
([IS]) taking 

m L = N(N- l)/2. (19) 
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IV. RESULTS AND DISCUSSION 

In what follows, we display the results obtained from 
exact diagonalization described in the previous section. 

In order to study under which conditions the Laughlin 
state is the ground state of the system a phase diagram 
is presented in Fig. [5J The fidelity between the Laughlin 
state and the GS of the system is plotted versus the laser 
intensity Uq /Er and the strength of the contact interac- 
tion g. For each value of Uq/ En, we have taken the max- 
imum possible rotation frequency, defined by Eq. (fT5"j). 
since it corresponds to the best condition to achieve the 
Laughlin wave function. The plot shows two separated 
phases. In one of them, the Laughlin state is perfectly 
obtained ((GS\^l) — 1)- It corresponds to high values of 
the laser intensity and strength of the interaction. On the 
contrary, for low values of Uo/Er and g the fidelity be- 
tween the GS and the Laughlin state is zero. This phase 
corresponds to other states with less angular momenta 
and not so correlated. 

Notice that the transition between the Laughlin and 
non-Laughlin phases is very abrupt. This is because 
the ground state in the non-Laughlin phase has an an- 
gular momentum smaller than the angular momentum 
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FIG. 2: Density plot of the fidelity between the GS of the 
system and the Laughlin wave function versus the laser in- 
tensity Uo and the strength of the contact interaction g for 
different number of particles (./V=3,4 and 5). In order to re- 
alize the diagram, the maximum possible rotation frequency 
Ql has been considered for each value of Uo- The dashed line 
represents the dependence of g on the confinement Uo for the 
Rubidium atoms case. Although the Uo needed to get the 
Laughlin as the ground state of the system increases enor- 
mously with the number of particles, we realize that for small 
values of N, the laser intensity required is achievable. 



G 



that Laughlin state requires. States with different angu- 
iar momentum are orthogonai and therefore the fidelity 
between them is zero. Nevertheless, when interaction is 
high enough and it is sufficiently favorable for the sys- 
tem to have total angular momentum L = N(N — 1), 
the Laughlin wave function is instantly achieved. This 
argument is illustrated in Fig. [3] for the case of N = 4 
particles. We can see how the system increases its angu- 
lar momentum depending on the confinement Uo and the 
interaction g. 

In Fig. [2j the shape of the boundary between the 
Laughlin and non-Laughlin phases can be easily ex- 
plained. According to the argument presented in the 
previous section, the configuration of the ground state 
of the system depends on the ratio between the strength 
of the interaction g and the energy difference of the single 
particle spectrum AE = E mL — Eq = ^m^im^ + 1) (see 
Fig. [IJ. In particular, we expect two regimes: 

• If g ^> AE, the GS is the Laughlin state, 

• while if g <C AE, then GS is a product state with 
all the particles with angular momentum 0. 

The transition, then, takes place between these two 
regimes, that is, when g ~ AE. This condition implies 
that 



V Uq 



(20) 



where f(N) ~ ra^ijn^ + 1) is a function that depends on 
N. ' 

In Fig. SI we plot the numerical data (<J jjf,g) corre- 
sponding to the points of the border between the Laugh- 
lin and non-Laughlin phases of Fig. [5] for N — 4 in order 
to see if Eq. (|2"0)> is fulfilled. We observe a perfect fit be- 
tween numerical data and Eq. (|20p and, performing a lin- 
ear regression, we determine the slope /(4) = 33.876(1). 
The same perfect agreement has been found for N = 3 
and N = 5 cases, with slopes /(3) = 10.779(1) and 
/(5) = 77.427(3) respectively. 

According to Eq. ([20]). we would expect that function 
f(N) scales as - N(N - 1) (N(N - 1) + 2). Neverthe- 
less, although /(3), /(4) and /(5) seem to scale in this 
way, three points are not enough to confirm this behavior. 

The dashed line in Fig. [5] expresses the dependence of 
g on the confinement Uo for the Rubidium atoms case 
(see Eq. ([2"Tj) L It shows the natural path that the system 
would follow increasing the intensity of the laser without 
improving the strength of the interaction by means of 
Feshbach resonance techniques. This illustrates the ad- 
vantages of confining ultra-cold gases in an optical lattice. 
The high confinement increases the scattering length of 
the interaction between the atoms, 
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FIG. 3: Density plot of the angular momentum of the GS of 
the system respect to the laser intensity Uo and the strength 
of the contact interaction g for iV=4. In order to realize the 
diagram, the maximum possible rotation frequency Ql has 
been considered for each value of Uo- 
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FIG. 4: Plot of the points (w j^,g) of the border between the 
two phases Laughlin and non-Laughlin of Fig. [5] for N = 4. 
We realize that Eq. [20] describes perfectly its behavior. 



which, at the end of the day, is responsible for the high 
correlation of the GS. 

Moreover, this enhanced interaction produces a larger 
energy gap which makes the optical lattice proposal more 
robust compared to the harmonic trap setup. This en- 
hancement in the energy gap of the excitation spectrum 
can alleviate some of the challenges for experimental re- 
alization of the quantum Hall state for ultra-cold atoms. 

In Fig. [5J we also realize that the larger the number of 
atoms of the system is, higher the laser intensity needed 
to achieve the Laughlin state. The reason for this is 
that for larger N, the maximum angular momenta al- 
lowed before the atoms escape from the trap is larger, 
mi = N(N — l)/2. This implies a smaller maximum 
rotation frequency £Il (see Eq. (|15[) L and the require- 
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merit of increasing the laser intensity Uq to compensate 
this effect. Let us note that increasing Uo we achieve 
both larger maximum rotation frequency and a stronger 
interaction g. 

It is interesting to point out that driving the system 
into the Laughlin wave function is not possible by sim- 
ply increasing the rotation frequency, since the symmet- 
ric shape of the potential conserves the angular momenta 
and the system would rest in a trivial non-entangled state 
with angular momentum zero. To avoid this, we should 
introduce some deformation to the wells of the lattice in 
order to break the spherical symmetry. For our optical 
lattice setup this can be achieved by introducing a cou- 
ple of electro-optical fast phase modulators whose effect 
in the rotating frame would be a new trapping potential 
with the form V p oc (w + Auj) 2 x 2 + u 2 y 2 . This modifica- 
tion would change our original Hamiltonian in a new one 
H + H e , where 

H < = I £ ^ m+2 b m + (m + l)blb m + h.c, (22) 

m 

where j3 m = y/ (m + 2)(m + 1) and e = Alo/uj is a small 
parameter. 

The previous perturbation (|22p leads to quadrupole 
excitations, so that states whose total angular momen- 
tum differ by two units are coupled, and the system can 
increase its angular momentum. In order to drive the 
system into the Laughlin state, we could follow some adi- 
abatic paths in the parameter space (f2, e) such that the 
gap along this paths is as large as possible, and in this 
way, keep the system in its ground state according to the 
adiabatic theorem [la ]. 

Finally, let us briefly discuss how to check that we have 
driven the system to the Laughlin state. First, notice 
that as the system is in a Mott phase, any measurement 
signal will be enhanced by a factor equal to the number 
of occupied lattice sites. 

According to Ref. [16], for any LLL state of a trapped, 
rotating, interacting Bose gas, the particle distribution 
in coordinate space in a time of flight experiment is re- 
lated to that in the trap at the time it is turned off by 
a simple rescaling and rotation. Thus, it is possible to 
measure the density profile of our state by accomplishing 
a time of flight experiment. Notice that by means the 
density profile it is possible to estimate the total angular 
momentum of the system. 

Although the angular momentum of the system can 
be measured, there are many states with the same angu- 
lar momentum that are not the Laughlin state. Thus, 
in order to really distinguish if the system is in the 
Laughlin state, the measurement of other properties is 



required. One possibility is the measurement of correla- 
tions. In Ref. [l3j, an interesting technique is proposed 
to measure the correlation functions g\ — (ifi (r)ijj(r')} 
and #2 = (■0^( r )V^( r ')' ! /'( r )' ! / ; ( r '))- These functions are 
very particular for the Laughlin state. For instance, 
52 oc \r — r'| 4 , since in the Laughlin state all the parti- 
cles have a relative angular momentum m r = 2 between 
them. Let us just mention the scheme of this technique. 
It lies in using a gas of two species that can be coupled 
via Raman transitions (hyperfine levels). First, all the 
atoms are in the same state and the system is driven to 
the Laughlin state. Then, an equal superposition of the 
two species is created by means of a 7r/2-pulse with the 
laser. Next, the atoms are shifted a small distance com- 
pared to the lattice spacing by moving the lattice of one 
of the species (as was proposed in Ref. [I?} and experi- 
mentally realized in Ref. [18:]). Finally, another 7r/2-pulse 
is performed in order to put all the atoms in the same 
state and the time of flight measurement is accomplished. 



V. CONCLUSIONS 

The method for achieving the Laughlin state in an op- 
tical lattice presented in Ref. [l3| has been studied. It 
has been shown that it is essential to take into account 
the anharmonic corrections. A quartic anharmonic term 
introduces a maximum rotation frequency that is smaller 
than in the harmonic case and, therefore, if this was not 
considered all the particles would be expelled from their 
single lattice sites. 

Although this more restrictive centrifugal limit makes 
more difficult to drive the system into the Laughlin state, 
since it requires higher values of Uq, for systems with a 
small number of particles, the Laughlin state is achiev- 
able experimentally. 

The shape of the boundary between the Laughlin and 
non-Laughlin phases in the phase diagram of the system 
respect to the confinement Uq and the strength of the in- 
teraction g plotted in Fig.[5]has been physically explained 
and analytically described. 



Acknowledgments 

The main part of this work has been realized under the 
supervision of B. Paredes. Let us also thank S. Trotzky, 
S. Iblisdir and M. Moreno for fruitful discussions. Finan- 
cial support from QAP (EU), MICINN (Spain), FI pro- 
gram and Grup Consolidat (Generalitat de Catalunya), 
and QOIT Consolider-Ingenio 2010 is acknowledged. 



[1] D. C. Tsui, H. L. Stormer, and A. C. Gossard, Phys. Rev. [2] R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

Lett. 48, 1559 (1982). [3] D. Yoshioka, The quantum Hall Effect ( Springer- Verlag, 



8 



Berlin, 2002). 

[4] X. G. Wen, Quantum Field Theory and Many-Body Sys- 
tems (Oxford University Press, Oxford, 2004). 
[5] A. Kitaev, Annals of Ph ysics 321, 2 (2006), URL 

|arXiv : cond-mat/0506438] 
[6] B. Dougot, M. V. Feigel'man, L. B. Ioffe, and A. S. 

Ioselevich, Phys. Rev. B 71, 024505 (2005), URL 

|arXiv : cond- mat/0403712l 
[7] L.-M. Duan, E. Demler, and M. D. Lukin, Phys. Rev. 

Lett. 91, 090402 (2003), URL |arXiv : cond-mat/021 0564 
[8] A. Micheli, G. K . Brennen, and P. Zoller, N at. Phys. 2, 

341 (2006), URL [arXiv : quant-ph/0512222 1 
[9] B. Paredes, P. Fedichev, J. I. Cirac, and P. Zoller, 

Phys. Rev. Lett. 87 , 010402 (2001), URL 

|arXiv: cond-mat/010325l1 
[10] B. Paredes, P. Zolle r, and J. I. Cirac, Phys. Re v. A 66, 

033609 (2002), URL |arXiv: cond-mat/0203061| 
[11] M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, 

A. S. De, and U. Sen, Adv. in Phys. 56, 243 (2007), 



URL arXiv:cond-mat/0606771 

[12] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008), URL IarXiv:0704.301ll 
[13] M. Popp, B. Paredes, and J. I. Cirac, Phys. Rev. A 70, 

053612 (2004), URL |arXiv : cond-mat/0405195 1 
[14] C. Becker, P. Soltan-Panahi, J. Kronjager, 

S. Dorscher, K. Bongs, and K. Sengstock (2009), 

URL larXiv: 091273 646 
[15] M. Born and V. A. Fock, Zeitschrift fur Physik a Hadrons 

and Nuclei 51, 165 (1928). 
[16] N. Read an d N. R. Cooper, Phys. Rev . A 68, 035601 

(2003), URL [arXiv : cond-mat/0306378 1 
[17] D. Jaksch, J. I. Cirac, P. Zoller, S. L. Rolston, R. Cote, 

and M. D. Lukin, Phys. Re v. Lett. 85, 2208 (2000), URL 

|arXiv : quant- ph/0004038 
[18] O. Mandel, M. Greiner, A. Widera, T. Rom, T. W. 

Hans en, and I. Bloch, Phys. Rev. Lett. 91, 010407 (2003), 

URL [arXiv : cond-mat/0301 169] 



